Radiative association and inverse predissociation 

of oxygen atoms 

J. F. Babb and A. Dalgarno 

Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138 

^" ! (February 5, 2008) 

OS 



> 

o 
in 



X 



Abstract 



The formation of O2 by radiative association and by inverse predissociation 



J> \ of ground state oxygen atoms is studied using quantum-mechanical meth- 

o 

ods. Cross sections, emission spectra, and rate coefficients are presented and 
compared with prior experimental and theoretical results. At temperatures 
0^ I below 1000 K radiative association occurs by approach along the 1 IT n state 

Q_i! of O2 and above 1000 K inverse predissociation through the B 3 £~ state is the 

dominant mechanism. This conclusion is supported by a quantitative com- 

r-j 

(\ parison between the calculations and data obtained from hot oxygen plasma 



spectroscopy. 
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I. INTRODUCTION 

Of the 81 molecular electronic states that two 0( 3 P) atoms can form, the repulsive 1 3 U U 
state can make an allowed electric dipole transition to the X 3 £~ ground state. Thus, the 
formation of an oxygen molecule through the direct radiative association reaction 

o( 3 p) + o( 3 p) - o 2 (i 3 n u ) - o 2 (x 3 s 9 -) + hv (i) 

can occur and it might be an important molecular oxygen creation process at low tempera- 
tures. In Fig. |T| the potential energy curves for the 1 3 IT U and X 3 E~ states of 2 are depicted. 
Wraight fl|] had suggested that the process (|J) might be responsible for the emission in the 
near-infrared of the nightglow of the planet. 

At higher temperatures another process — inverse predissociation, an indirect mechanism 
of association — can occur with the excited B 3 E~ state serving as an intermediate step in 
the formation of an oxygen molecule 

0( 3 P) + 0( 3 P) - O* - 2 (B 3 E~) - 2 (X 3 S 9 ") + hv, (2) 

where 0* 2 denotes any one of the four states 1 1 I1 U , 1 3 H U , 1 5 IT U , or 2 3 £+ that dissociate to 
0( 3 P) + 0( 3 P) and cross the B 3 S~ state. The predissociative coupling is due to spin-orbit 
(any of the four states) or electronic-rotational interaction (1 3 U U state) ||. Reaction (0) may 
be characterized by an activation energy that is the energy with respect to the 0( 3 P) + 0( 3 P) 
asymptote at which the particular state 2 crosses the B 3 S~ state. The potential energy 
curve of the B 3 S~ state, which dissociates to an 0( 3 P) atom and an 0( X D) atom is shown 
in Fig. [I]. Because the strength of the electronic-rotational coupling increases with the 
nuclear rotational angular momentum, while the strength of the spin-orbit coupling is nearly 
independent of it || , the 1 3 U U state is expected to be the most important channel for 
inverse predissociation. In a shock tube study of the emission from recombining oxygen 
atoms in a non-equilibrium shocked ozone-argon gas mixture at temperatures between 2500 
and 3800 K Myers and Bartle obtained absolute intensity measurements of the spectrum, 



and attributed the emission to the formation of oxygen molecules from ground state oxygen 
atoms through the inverse predissociation process 

0( 3 P) + 0( 3 P) - 2 (1 3 n u ) - 2 (B 3 S,,;) - 2 (X 3 S;) + hv. (3) 

The emission has been observed in several later experiments. It was seen in a silent electric 
discharge in 2 by Weisbeck and Volkner || and Sharma and Wray || measured emission 
intensities at 230 and 325 nm for temperatures in the range 2800-5300 K in shocked oxygen 
and oxygen-noble gas mixtures. Sharma and Wray || confirmed that the excited levels of 2 
were populated by inverse predissociation but their emission spectra appear to be modified 
by vibrational redistribution. Similar conclusions were reached by Wray and Fried |7j who 
measured the emission intensities in an arc jet in oxygen and and oxygen-like gas mixtures. 
Hoffmann and Neiger 0,[| measured the absolute intensity of the emission from an equilib- 
rium oxygen plasma at 8800, 9140, and 9420 K, finding a peak at about 260 nm which they 
attributed to inverse predissociation (|3]). We return to a discussion of prior work in Sec. |V|. 
In this paper, we study quantum-mechanically the direct radiative association process (|XJ) 
and the inverse predissociation process (^). Cross sections, emission spectra, and rate coeffi- 
cients are determined for temperatures from 15 to 15 000 K. Improvements in the knowledge 
of the 1 3 U U and B 3 £~ state potentials make it possible to get qualitative agreement with 
many experiments for the emission spectra and quantitative agreement for rate coefficients 
measured by Myers and Bartle ||. 

II. MOLECULAR STATES 



In Ref. |10 we constructed two 1 3 I1 U potential energy curves that were identical except 
in the range of internuclear distances 2.4 < R < 3.1 (we shall use atomic units in this section 
and Sec. [TTI| ). The first curve was designed to be consistent with the level shifts arising from 



predissociation of the B E u state, as analyzed in Ref. [11] where an internuclear distance 



for the B 3 S U - 1 3 U U crossing R x = 2.700 was obtained, and which used ab initio energy 



calculations [12] for 3.1 < R < 10. The second curve we constructed simply used the ab 



initio calculations |12| for 2.5 < R < 10. A recent analysis pi of the photoabsorption cross 



section measurements [0, an analysis presumably more accurate than that of Ref. [jTT 



yielded a different crossing point of R x = 2.748, which is consistent with the crossing point 
of the cubic spline fit |T3[ to the calculations of \Y1\ at R = 2.5 and 2.8. We have therefore 
used the second curve we constructed in Ref. |10| here. As we have discussed |l(| the 



potential curve was supplemented for R < 2.4 by the adjusted ab initio energies of [|15|] and 
was constructed to have the proper asymptotic long-range behavior. 



The adopted transition dipole moment function D(R), defined in p6| , for the transition 
from the 1 3 fI M to the X 3 E~ state was the semi-empirical determination of |17] for 1.5118 < 
R < 2.5 and the ab initio calculation of [15| for 2.6 < R < 5. The data were smoothly 
connected using cubic splines. For R < 1.5118 the function D(R) = 0.1794 — 0.01945-R was 
used and for R > 5 the function D(R) = 1.2625/ R 3 . The X 3 £~ state potential energy curve 
was constructed using RKR data for 1.79480 < R < 4.0840 and fitted to the long-range form 
—17 .57 / R 6 for R > 5 as described in Ref. I£|] , but with a short-range form 



V(R) = 5757.0053 exp(-5.7925 J R). 



(4) 



The adopted X 3 E~ and 1 3 U U potential curves, which correlate with the separated atoms 
0( 3 P) + 0( 3 P), are shown in Fig. |l|. The B 3 S~ potential curve of |1| and the semi-empirical 
B 3 E~-X 3 £~ transition dipole moment of [0 as described in |1| were utilized. In Fig. |l|the 
adopted B 3 S^ potential is shown. It correlates to the separated atom limit 0( 3 P) + 0( 1 D, 
which lies 0.07229 au (1.967 eV) above 0( 3 P) + 0( 3 P). 



III. DIRECT RADIATIVE ASSOCIATION 



A fully quantum-mechanical theory can be formulated and used to calculate the cross 
sections at various energies and rate coefficients at various temperatures for direct radiative 
association ([I]). Let E be the energy of relative motion of the atoms approaching in the 



1 3 U U electronic state. The cross section for a transition to a bound vibration-rotation level 
E V N N i, of the X 3 £~ molecular ground state is given in atomic units by |20| . |2T | 



A 2 *3 

&v"n"n{E) = — — (E + E v ii N n) Sjn j\M N n N \ (5) 

9 fyE 

where [i = 14583.10 is the reduced mass for 16 C>2, g\ = 3 is the statistical weight of the 1 3 U U 

state (the nuclear spin is zero so only one of each lambda-doubling level is populated) and 

g — 81 is the statistical weight of the 0( 3 P)-0( 3 P) pair, a is the fine structure constant, 

Sjn y j is the Honl-London factor, and 

roc 

M N ,,, N = dR?P v „ N „(R)D(R)xN(R), (6) 

Jo 

with tfj v "Ni'(R) the normalized bound state wave function and xn{R) the energy-normalized 

wave function for the partial wave with angular momentum N of the continuum 1 3 II n state. 

The wave functions ^V'iV" and xn were calculated by integrating the Schrodinger equation 

using the Numerov method and xn was matched to the asymptotic form 

Xn(R) ~(2fi/nk) 1/2 sm(kR + 5 N -±NTT), R ~ oo, (7) 

with k = (2/i.E) 1 / 2 and 5^ the phase shift. Due to the presence of identical nuclei in the 
diatom and the absence of nuclear spin, the X 3 S^ state is allowed only odd values N" . 
Let J = |N + S| be the magnitude of the total angular momentum for the l 3 Il n state, 
with S the total spin angular momentum and |S| = 1. Assuming that Hund's case (b) 
coupling applies, we have for iV 3> 1, J ~ N (in an unrelated study, case (b) was found to 



be a satisfactory approximation to intermediate coupling [|II|) and then, summing over the 



rotational branches [22| for the allowed fine structure transitions for a given value of iV we 



obtain the line strength values given in Table |, which possess the sum rule property p3|j24 



J2Sn»,n = (2S + 1)(2N + 1). (8) 



N" 



In carrying out the computations iV" was required to be odd, but N could be odd or even. 
The total cross section, summed over allowed transitions between partial waves N and 
final states v"N", 



a ( E ^ (ozmU S °v>n»n{E) (9) 

\m -r i-) v n N „ N 

is presented in Fig. § at various energies. Some structure, due to shape resonances, is 
apparent in the cross section at low energies. The cross section drops off rapidly for E > 0.4 
due to the loss of Franck-Condon overlap at small R. In computing the total cross sections 
values of N up to 15 were sufficient for E < 2.0 x 10~ 5 , up to 30 for £<3.0x 10~ 5 , up to 
100 for E < 0.015, and for E > 0.015 values of iV up to 188 were used. Since the X 3 E~ state 
has 45 vibrational levels each total cross section at a given E was a sum of many partial 
cross sections. Nevertheless, the computations could be carried out readily. 

The cross sections obtained at various values of E were used to obtain the rate coefficient 
at various temperatures by averaging over a Maxwellian velocity distribution, 

1/2 

kT 

The values are given in column 2 of Table |H| The rate coefficient is approximated to within 
5% for T < 400 K by the expression 

a D (T) « 9 x i(r 27 T- a23 exp(T/70), T < 400 K. (11) 

Wraight [|IJ estimated semiclassically the rate coefficient for the direct process (]!]) and 
found 

a w (T) ~6x 10~ 24 T L58 cmV 1 , (12) 



^ = H^r) />(*>«*• HO) 



for 200 K < T < 2000 K. He advised caution in applying ( 120 to high temperatures. 
Expression (|i"2|) was adopted for astrophysical purposes in the tabulation of Prasad and 
Huntress p5| and in the UMIST data base |26| . The present computations show that Eq. (|12"D 



drastically overestimates the radiative association rate coefficient at low temperatures. 

IV. INVERSE PREDISSOCIATION 

For inverse predissociation, (|3[), the rate coefficient at temperature T is given by the 
Breit-Wigner expression [p^ -|29l 



^) = -(^f 2 E E i™* + 1) /; ( ^m ^f / ^ S «xp(- W^), (13) 

1 ; 2 V^ T / ^^ ; r r (*/,iv) + r d (t/,iV') PV vJV/ '' l ; 

where g-Q = 3 is the statistical weight of the B 3 E~ state, g = 81 is the statistical weight 
of the 0( 3 P)-0( 3 P) pair, E v > N i is the energy of the level v'N' of the B 3 E~ state measured 
with respect to the 0( 3 P)-0( 3 P) asymptote, T T (v',N') and T d (v',N') are, respectively, the 
radiative and predissociative widths of level v'N', with 

T T (v', N') = hA(v', N')=hJ2 A(v'N'v"N"), (14) 

v"N" 

and A(v'N'v"N") is the Einstein ^-coefficient. The sum in Eq. ( fL3|) is taken over levels v'N' 
with energy greater than the energy of the crossing point. Analysis |13| of the predissociation 
in the B 3 E~ state indicates that the 1 3 U U state crosses at an energy of 0.0435 relative to the 
0( 3 P)-0( 3 P) asymptote |IT| or about 13 000 K. For the B 3 £- state, typically T d {v',N') > 



r T (v',N') g, so that 

o / o_ \ 3 / 2 

a{T) W 81 (^rj ?^ (2iV' + l)r r K,iV')exp(-^ JV ,AT). (15) 

In evaluating (|T5|) we made the approximation 

A{v'NWN") « A(i/0/0)[(^ - E V „ N „)/{E V , Q - E vlll )\\ (16) 

with N' ~ iV" leading to the formula 

x J2( 2N> + ^[(Evn' - E V „ N ,)/{E V , Q - E v/ll )] 3 eM- E VN'/kT). (17) 

N> 

The intensity of emission from a transition from level (v', N') into level (v", N") is then 
given by the expression 

l v « /^^ J vv'iV''|[(27r//iA;T) 3 / 2 (2iV , + 1) 

x exp(-E„/ iV //A;T)M(t; / 0i; // 0)[( J E^ 7V / - E v ,, n ,)/{E v ,q - E v , fl )] 3 . (18) 

With E(v'N'v"N") = hv v >N>v»N» and JS^ both in cm" 1 , T in K, and A(v'0v"0) in s -1 , we 
can express I v in W-cm 3 as 

J„ = A(4.7 x 10~ 45 )£(t/NV'iV")(l/T 3/2 )(2iV' + l)A(uWO) exp(-^ 7V '/0.695T). (19) 



V. RESULTS AND COMPARISON WITH EXPERIMENT 



Using a program of Allison [31,19| and the potentials and transition dipole moments of 
Sec. [TT] the radiative widths were determined. They are in good agreement with other com- 
putations [31,32]. In Table [Tfl| we present the values of the emission transition probabilities 
A(v', 0) of the B 3 S~ state vibrational levels. In Fig. |3| the calculated emission spectra arising 
in inverse predissociation at 3000 K are shown in W-cm 3 /nm along with the measurements 
of Myers and Bartle 0] at 3030 K; the experimental points and error bars were interpolated 
from absolute emission-intensity spectral distribution curves. Our calculations support the 
suggestion that the measurements are in error for wavelengths A ;$ 300 nm, a possibility 
noted by Myers and Bartle Q and Sharma and Wray ||. The measurements of Sharma 
and Wray |J and Wray and Fried |7| appear to be affected by vibrational redistribution 
and their spectral data can be better interpreted in terms of an equilibrium model of the 
vibrational level populations of the B 3 E~ state |3^| . 

In the measurements of Hoffmann and Neiger near 9000 K a strong broad continuum 
feature was observed at around 260 nm. Our calculated emission spectra from inverse pre- 
dissociation are consistent with the experimental observations. Measurements at higher 
wavelength resolution should reveal structure arising from the Schumann-Runge bands 
populated by inverse predissociation. We find an absolute emission intensity of about 
5 x 10~ 39 W-cm 3 /nm at 250 nm and 9400 K. From the measured absorption coefficient 
of K = 2 x 10~ 5 cm -1 H, using Wien's radiation law P,Q, the length 10 cm of the exper- 
imental arc column, and estimating the number density |§ of O atoms of about 10 18 cm -3 
from Fig. 5 of Ref. @, we obtain an experimental value of 2.5 x 10 -39 W-cm 3 /nm at 250 nm 
and 9400 K. The agreement to within a factor of two is remarkably close. 

A rate coefficient 

«e xp (T) = 7.3 x 1(T 18 exp[(-14540 ± 1100) /T] cmV 1 , 280 < A < 500 nm, (20) 
for the total radiative combination of oxygen 



+ 0^0 2 + hu (21) 

between 2500 K and 3800 K was deduced by Myers and Bartle 0. In Table [TT] values of 
&Exp{T) are given for various temperatures. Our values for the rate coefficient for inverse 
predissociation a\p(T) are included in Table |TJ. The rate coefficient is fitted to within 10% 
for T < 5000 K by the expression 

aip(T) ~ 6.4 x 10~ 18 exp(-13 455/T) cmV 1 , T < 5000 K, (22) 

and to within 5% for 5 000 < T < 15 000 K by 
arp(T) ~ (-5.6 x 10~ 19 + 2.4 x 10" 22 T - 9.6 x 10~ 27 T 2 ) cmV 1 , 5 000 < T < 15 000 K. 

(23) 

In Fig. f| the calculated direct «d and inverse predissociation a IP rates are presented and 
compared with the «Exp- The agreement is very good, confirming that the formation of 
molecular oxygen from hot oxygen atoms at low densities proceeds by inverse predissociation 
through the B 3 S~ state following approach along the 1 3 U U state. 
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TABLES 
TABLE I. Line strengths for transitions from the l 3 II n to the X 3 S~ state. 

J" Branch Sj»,j/(2S + 1) 
J -I R J+l 

J Q 2J + 1 

J + l P J 
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TABLE II. Rate coefficients in cm 3 s , as a function of temperature T, in K, for direct radiative 
association process (j!]), inverse predissociation (||), and the experimentally derived value for the 
total process (|2l|). Numbers in square brackets represent powers of 10. 







Inverse 


Experiment, Ref. jl] 


T{K) 


Direct 


Predissociation 


a Exp , Eq. © 


15 


6.0[-27] 






30 


6.4[-27] 






50 


7.4[-27] 






100 


1.3[-26] 






200 


4.6[-26] 






300 


1.7[-25] 






400 


6.5[-25] 






500 


1.9[-24] 






600 


4.4[-24] 


1.1 [-27] 




700 


8.3[-24] 


2.7[-26] 




800 


1.4[-23] 


3.1[-25] 




900 


2.1 [-23] 


2.0[-24] 




1000 


2.9[-23] 


9.2[-24] 




1200 


4.7[-23] 


8.9[-23] 




1400 


6.9[-23] 


4.5[-22] 




1600 


9.2[-23] 


1.5[-21] 




1800 


1.2[-22] 


3.9[-21] 




2000 


1.4[-22] 


8.4[-21] 


5.1[-21] 


2200 


1.6[-22] 


1.5[-20] 


9.8[-21] 


2500 


2.0[-22] 


3.2[-20] 


2.2[-20] 


2700 


2.2[-22] 


4.8[-20] 


3.3[-20] 


3000 


2.5[-22] 


7.8[-20] 


5.7[-20] 
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4000 


3.7[-22] 


5000 


6.1[-22] 


L0000 


1.7[-20] 


L5000 


7.3 [-20] 



2.2[-19] 1.9[-19] 

3.9[-19] 4.0[-19] 

9.0[-19] 
9.3[-19] 
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TABLE III. Transition probabilities A(v',0) for emission from the B 3 S U to X 3 £ states in 
units of 10 7 s _1 for various values of the B 3 S~ vibrational quantum number v' with N' = 0. 



v' 

1 
2 
3 
4 
5 
6 
7 
8 
9 
10 



A(v',Q) 



v 



A(v',0) 



1.816 
1.916 
2.006 
2.078 
2.130 
2.157 
2.157 
2.129 
2.073 
1.992 
1.883 
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12 
13 
14 
15 
16 
17 
18 
19 
20 
21 



1.750 

1.586 

1.393 

1.184 

0.9709 

0.7833 

0.6162 

0.4710 

0.3414 

0.2254 

0.1195 
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FIG. 1. Adopted 1 3 II U , B 3 S", and X 3 S" potentials. 
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FIG. 2. Cross section a{E) @ in au (oq) for the direct radiative association process (|T|) as a 
function of energy E in au (e 2 /ao). 
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FIG. 3. Absolute emission intensity at various wavelengths calculated for inverse predissociation 
at 3000 K (solid line), compared to experimental measurements of Myers and Bartle at 3030 K 
(points with error bars connected by the dotted line). The calculated spectrum has been averaged 
over bins of width 5 nm. The experimental data are average values over a spectral interval the 
width of which is represented by the horizontal span of the error bar. 
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FIG. 4. Comparison of the calculated direct radiative association (solid line) and inverse pre- 
dissociation (dotted line) rate coefficients and the experimental rate coefficient (points with error 
bars) of Myers and Bartle. 
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